Mon. Not. R. Astron. Soc. 000, [JHTS] (2002) Printed 22 February 2013 (MN WI^ style file v2.2) 



Kinematics and chemistry of the hot core in G20.08-0.14N 



Jin-Long Xu^'^* and Jun-Jie Wang^ 

^National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China 
" NAOC-TU Joint Center for Astrophysics, Lhasa 850000, China 



4,2 



o 

(N 
(N 

^' 

o 



> 

O 

cn 

o 

m 



22 February 2013 



ABSTRACT 

We present Submillimeter Array observations of the massive star-forming region 
G20.08-0.14N at 335 and 345 GHz. With the SMA data, 41 molecular transitions 
were detected related to 11 molecular species and their isotopologues, including SO2, 
SO, C^^S, NS, C^^O, SiO, CH3OH, HC3N, H13C0+, HCOOCH3 and NH2CHO. In 
G20.08-0.14N, 10 transition lines of the detected 41 transition lines belong to SO2, 
which dominates the appearance of the submillimeter-wave spectrum. To obtain the 
spatial kinematic distribution of molecules in G20.08-0.14N, we chose the strongest 
and unblended lines for the channel maps. The channel maps of C'^'^S and SiO, to- 
gether with their position- velocity diagrams, present that there are two accretion flows 
in G20.08-0.14N. Additionally, SiO emission shows a collimated outflow at the NE-SW 
direction. The direction of the outflow is for the first time revealed. The rotational 
temperature and the column density of CH3OH are 105 K and 3.1 x 10^^ cm~^, re- 
spectively. Our results confirm that a hot core is associated with G20.08-0.14N. The 
hot core is heated by a protostar radiation at it center, not by the external excitation 
from shocks. The images of the spatial distribution of different species have shown that 
the different molecules are located at the different positions of the hot core. Through 
comparing the spatial distributions and abundances of the molecules, we discuss pos- 
sible chemical processes for producing the complex sulfur-bearing, nitrogen-bearing 
and oxygen-bearing molecules in G20.08-0.14N. 
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1 INTRODUCTION 

Massive stars are formed in dense molecular clouds. Massive 
star formation also has a significant effect on the chemistry 
of th e surrounding molecular clouds (Ivan Dishoeck fc Blakd 
Il998h . Hot cores are the formation sites of massive stars. 
The hot cores are defined as compact ( 5C 0.1 pc, n ^ 10^ 
cm~'^), relatively hi gh temperature (Tk ^ 100 K) cloud cores 
jKurtz et al. I 2000f). whose phase is thought to last about 
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van Dishoeck &: Blakd 1 19981 1. When the hot core is 



formed, the central massive protostar can produce the ioniz- 
ing radiation and the associated outfiow can produce shocks, 
hence the hot cores represent the most chemically rich 
phase of the ma ssive star formation often associated with 
UC HII regions JCesaroni et al."l992': Hatche U et al.|[l99i : 
iGarav fc Lizand Il999l : IChurchwell 2002). The high abun- 
dances of the organic molecules in hot cores are consequently 
attributed to grain-surface chem istry and mantle eyapora - 
tion processes (|van Dishoeck fc Blake 199^ : lLiu et al.ll200ll ). 
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Because of the compact and dense nature of hot cores, single- 
dish observations with large beam sizes are not sufficient to 
explore the dense cores and detailed kinematics. The in- 
terferometer observations at submillimeter wave can filter 
out the extended diffuse components, then the detailed dy- 
namical processes and chemical conditions of hot cores can 
be revealed. A number of line observations at submillime- 
ter/millimeter wavelengths were previously used to explore 
the molecular composition of hot c ores ( Bcutho r ct al. 20051 : 
iGoddi et al.l l2009l : lOin et~all I2OI0I ) , but more sources need 
to be added to the inventory of studied hot cores before their 
chemical evolution can be understood. 

G20.08-0.14N is a massive star-form ing region. It is 
approxim ately at a d istance of 12.3 kpc (|Fish et al.ll2003l : 
[Anderson et al.l |2009| ) , corres ponding to a bolometric lu - 
minosity of about G.GxlO'^LM llGalvan-Madrid et all bOOgl ). 
Previous radio continuum observations at centimeter wave- 
lengths suggest that G20.08-0.14N has t hree UC and 
HC HII regions (|Wood fc Churchwelll [l989l ). In addition, 
single-dish observations of molecular li nes show the signa- 
tures of infall, accretion, and outfiow (jKlaassen fc WilsonI 
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12007 20081). Moreov er. H2O (iHofner fc C hurchwelll Il996l ). 
OH (Iho et allll983l'). CH3O H (|Walsh et al..,1998. ') and NH3 
iGalvan-Madrid et al.l [20091 ') masers in G20.08-0.14N were 
revealed in some observations. These signatures indicate ac- 
tive massive star formation in this region. The observed 
CH3CN transitions indicate t hat a hot molecular core i s 
associated with G20.08-0.14N (|Galvan-Madrid et al.ir2009l ). 
Thus, G20.08-0.14N provides us with an opportunity to 
study the physical and chemical conditions of massive star- 
forming processes. 

We have carried out multiline observations toward the 
massive star- forming region G20.08-0.14N with the Submil- 
limeter Array. Various molecular lines are used to investi- 
gate the physical and chemical processes of G20.08-0.14N. 
In Section 2, we summarize the observations and data re- 
duction. In Section 3, we give the general results. In Section 
4, we present the data analysis for deriving rotation temper- 
atures, column densities, and abundances of various species 
relative to H2 and implications for chemistry. In Section 5, 
we summarize our main conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations toward G20.08-0.14N were carried out with 
the SMA on 2009 May 20, at 335 (lower sideband) and 345 
GHz (upper sideband). The data are from SMA archivsQ- 
The two sidebands of the SMA covered frequency ranges 
of 335.6-337.6 GHz and 345.6-347.6 GHz, respectively. The 
total observing time is 9.38 hours. The phase track cen- 
ter was R.A.(J2000.0)= 18^28^10=. 30 and Dec.(J2000.0) = 
— 11°28'47".8. The typical system temperature was 186 K. 
The spectral resolution is 0.812 MHz, corresponding to a 
velocity resolution of 0.7 km s~^. The bright quasar 3C273 
was used for bandpass calibration, while absolute flux den- 
sity scales were determined from observations of Callisto (15 
Jy). QSO 1733-130 and QSO 1751+096 were observed for 
antenna gain corrections. The calibration and imaging were 
performed in Miriad. The continuum image was constructed 
from the line- free channels. The spectral cubes were con- 
structed using the continuum-subtracted spectral channels. 
Self-calibration was performed on the continuum data. The 
gain solutions from the continuum were applied to the line 
data. The synthesized beam size of the continuum was ap- 
proximately 2".02 X 1".15 with a P.A. =72°. 0. 



3 RESULTS 

3.1 Continuum Emission at 0.9 mm 

Figure 1 shows the 0.9 mm continuum map of G20.08-0.14N 
obtained by the SMA observations. Galvan-Madrid et al. 
(2009) resolved the G20.08-0.14N system into three compo- 
nents with the VLA at 1.3 cm. Each component represents 
an H II region marked in our Figure 1 (,Wood fc Churchw ell 
Il989l ). Particularly, the HII region A is the brightest and 
closest to the peaks at 1.3 mm and 0.9 mm. By using a 
two-dimensional Gaussian fit for the continuum emission, we 
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Figure 1. Continuum map toward G20.08-0.14N at 0.9 mm. The 
contours are at -4, 4, 8, 12, 20, 28, 36, 44 and 49 a. The rms 
noise level is 0.03 Jy beam~^(l(T). The synthesized beam (2.0" X 
1.2") with P.A.= 72° is shown in the lower right corner. " * " 
indicates the position of three HII region (Wood & Churchwell 
1989; Galvan-Madrid et al. 2009). H2O maser is shown with blue 
filled triangle. The velocity component of the blue contours is 
from 32.5 km s~^ to 36.0 km s~^, while the velocity component 
of the red contours is from 47.0 km s~^ to 51.7 km s~^, whose 
levels are the 20%, 40%, 60%, 80%, and 100% of the peak values. 



obtained that the total flux density is 2.72±0.03 Jy, the de- 
convolved source size is 1.83" x 1.22" (P.A.=72°), the peak 
position is R.A.(J2000)=18''28™10."307 (AR.A.=±0.01"), 
decl.(J2000)=-ll°28'47."846 (Adecl.=±0.0l") with an in- 
tensity of 1.33 ± 0.02 Jy beam~^. The peak position of the 
0.9 mm continuum emission is coincident with HII region 
A within the unc ertainty. An H2O maser has been detected 
in G20.08-0.14N jHofner fc Churchwelll Il996h . In Figure 1, 
the H2O maser is shown with the black filled triangle. The 
maser is associated with 0.9 mm continuum emission, but 
offset from its the peak position. 



3.2 Molecular Line Emission 

Molecular lines were identified following spectral-line cat- 
alogs: (1) Cologne Databa se for Molecular Spectroscopy 
ijCDMS. Miiller et all l2005l ). (2) Mole cular Spectroscopy 



database of Jet Propulsion Laboratory (|JPL. Pickett et al 
1998"), and (3) SPLATALOGUE line catalogj^l jRemiian 
2007,). 



The obtained SMA 4 GHz spectrum is shown in Figure 2 
and a full list of identified lines presented in Table 1. We de- 
tected a total of 41 transitions, which include 35 transitions 
from 11 species and their isotopologues, and 6 unidentified 
transitions marked with U. The identified species contain 
simple linear molecules as well as complex oxygen-bearing, 
nitrogen-bearing, and sulfur-bearing molecules. Many of the 
strongest lines in the spectrum are from diatomic molecules 



|http://www.cfa.harvard.edu/rtdc/data/search.html| 



|http: / /www.splatalogue.net | 
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Figure 2. SMA spectra extracted from the line data-cubes in the image domain toward the G20.08-0.14N. The spectra are averaged 
over one beam centered on the peak of the continuum. U marks the unidentified Hues. 



of the most abundant elements. The ^ CO J=3-2 line at 
345.938 GHz is the strongest line in our passband, but 
its spectral profile displays absorption features as well. Be- 
cause the CO absorption features are at the sam e veloc - 
ities as the HI absorption features of iFish et al.l l|2003h . 
iGalvan-Madrid et al.l (|2009l ') explained that these features 
may be caused by foreground gas that is not associated 
with G20.08-0.14N but rather to intervening Galactic spi- 
ral arms. In addition, sulfur- bearing molecules such as SO2 
and C'^'^S display also strong emission lines. The line param- 
eters are presented in Table 1. The first column is the name 
of molecular species. The second and third columns list the 
transition and rest frequency of the molecules. The fourth 
column lists the upper level energy of each transition, and 
the fifth column is the product of the fine strength and the 
square of the relevant dipole moment. The sixth, seventh, 
and eighth columns list the central line velocity, the peak 
intensity, and full width at half-maximum (FWHM) derived 
from Gaussian fitting to the line profiles, respectively. The 1 
a integrated-intensity noise level of each transition is given 
in the ninth column. 

To obtain the spatial kinematic distribution of 
molecules in G20.08-0.14N, we chose the strongest and un- 
blended lines for the channel maps if multiple lines for one 
species were detected. Here we only show the channel maps 
of C'^*S and SiO in Figures 3-4. The contours are made at 
the velocity channels with interval 1 km s^^, which begin at 



4 a. In Figures 3-4, we present the full velocity range of each 
molecular emission. The three HII regions are shown with 
"*" symbols. Since the HII region A is coincident with the 
peak position of the 0.9 mm continuum (within the uncer- 
tainty), the symbol of the HII region A can also represent the 
peak position of the 0.9 mm continuum. The sy stemic veloc- 
ity iVsYs) for G20.08-0.14N is abo ut 42 km s"^ (|Plume et all 
ll992l : lGalvan-Madrid et al.ll2009h . From Figures 3-4, we can 
see that both molecules show a velocity gradient across the 
0.9 mm continuum along the northeast-southwest direction. 
To further determine the peak move of the 0.9 mm con- 
tinuum, we made a dashed line across the HII region A of 
four panels in Figures 3-4. Also, the peak position of the 
continuum above the dashed line is redshifted, while the 
blueshifted emission is under the dashed line, further con- 
firmed that there is a velocity gradient in G20.08-0.14N. 
According to the channel maps of each species, we made 
the integrated intensity map of each species (Figure 5). In 
Figure 5, the emission peaks of CH3OH and HCOOCH3 are 
associated with the continuum peak. The emission peaks of 
"CH3OH and NH2CHO are located to the northeast of the 
continuum peak, while the emission peaks of HC3N, C^^O, 
C^^'S, NS, and SiO are located to the southwest of contin- 
uum peak. Three sulfur-bearing molecules (SO2, '^''802, and 
SO) are situated to the northwest of continuum peak. Addi- 
tionally, the emission of H^'^CO^ shows two molecular cores 
around the HII region A. 
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Table 1. Observed Parameters of Each Line 



Molecule 


Transition 


Rest Frequency 


Eu 




Ay Integrated-intensity rms Notes 






(MHz) 


(K) 


(debye^) (km s ^) (Jy beam 


^) (km s"-"^) (Jy beam"-"^ km s^^) 



SO2 


233 21 22 


336089.23 


276.0 


28.41 


41.8(0.1) 


4.5 


.0.1) 


8.8(0.1) 








-^^7, 9"-^ "'^6, 12 


336669.58 


245.1 


4.34 


41.7(0.1) 


3.5 


0.1) 


8.9(0.2) 








-^^4,12"-'^^3,13 


346523.86 


164.5 


23.10 


38.5(0.1) 


7.8 


0.1) 


13.5(0.2) 




1 




"^^1 , 19 ~"^^0 , IS 


346652.16 


168.1 


41.98 


41.6(0.1) 


7.0 


0.1) 


10.1(0.1) 


8.0 






2^1,19~-^'^2,18^2 — ^ 


336760.70 


962.1 


21.72 


42.2(0.2) 


1.2 


'0.1) 


7.2(0.4) 




2 




"^22|ll~Hl^l0^2 — ^ 


337191.5 


_ 


_ 


35.9(0.1) 


2.2 


'0.1) 


7.9(0.2) 




3 




19i ig-18i_ig?;2=l 


346379.19 


930.6 


43.22 


42.0(0.2) 


1.6 


'0.1) 


9.4(0.3) 








I84 ^4-183 i5t'2=l 


346591.78 


960.8 


27.65 


41.9(0.1) 


1.9 


'0.1) 


9.2(0.3) 




4 


CH3OH 




336438.22 


488.2 


2.38 


42.3(0.3) 


0.6 


'0.1) 


4.9(0.5) 








12i_44A--12o 12^+ 


336865.15 


197.1 


22.88 


42.6(0.1) 


4.8 


'0.1) 


6.5(0.1) 


3.2 








337135.86 


61.6 


0.25 


42.2(0.2) 


1.2 


'0.1) 


6.7(0.4) 








-'-^1,15^ "-^-^2,14-^ 


345903.97 


332.5 


7.13 


42.1(0.1) 


2.4 


'0.1) 


6.0(0.2) 










346202.77 


115.2 


0.50 


41.6(0.1) 


2.3 


'0.1) 


6.5(0.2) 










337297.44 


390.0 


5.55 


41.8(0.1) 


1.7 


'0.1) 


6.4(0.3) 








73 5A-63 4Aj;j=1 
7= jA+-6'nA+Vf=l 

0,0 , jJ t 


337519.1 


323.9 


4.6 


36.2(0.8) 


1.0 


0.2) 


6.4(1.8) 








337546.05 


485.4 


2.74 


42.1(0.2) 


1.0 


0.1) 


4.2(0.3) 




5 


HC3N 


37-36 


336520.08 


306.9 


55.25 


41.6(0.1) 


4.2 


0.1) 


8.4(0.1) 


2.0 






38-37 


345609.02 


323.5 


529.13 


41.7(0.1) 


3.8 


0.1) 


8.5(0.2) 








37-36e 


337335.35 


1025.3 


513.34 


35.4(0.2) 


2.4 


'0.1) 


9.9(0.4) 




6 




38-37e 


346455.73 


645.1 


526.17 


42.7(0.2) 


2.2 


0.1) 


11.3(0.4) 








38-37/ 


346949.12 


645.6 


526.19 


40.3(0.2) 


2.4 


0.1) 


15.6(0.4) 




7 


HCOOCH3 


2'''g_ig^-26g jgA 


336028.15 


277.4 


64.03 


39.7(0.5) 


0.4 


'0.1) 


7.4(1.0) 








275,22^-265 21 


347478.24 


247.3 


69.17 


37.3(0.2) 


1.2 


'0.1) 


9.1(0.5) 


1.1 






^4,2-^3,3 


345651.30 


51.8 


4.49 


42.5(0.1) 


2.5 


'0.1) 


7.7(0.2) 








17 -17 

-^'4,14 -^'3,15 


345929.29 


178.8 


24.49 


43.1(0.1) 


2.6 


'0.1) 


11.0(0.3) 


2.0 




NH2CHO 


-'^^2,15--'^^2,14 


336136.88 


149.7 


617.20 


43.5(0.2) 


1.0 


'0.1) 


4.8(0.4) 


0.5 




SO 


lOii-lOio 


336553.75 


142.8 


0.29 


41.8(0.1) 


2.7 


'0.1) 


8.2(0.2) 


2.1 




Ci^O 


3-2 


337061.10 


32.4 


0.04 


41.7(0.1) 


2.7 


'0.1) 


12.1(0.2) 


5.2 






7-6 


337396.46 


50.2 


25.57 


41.8(0.1) 


5.1 


.0.1) 


7.2(0.1) 


5.2 




CO 


3-2 


345795.99 


33.2 


0.04 












8 


NS 


15/2, 17/2 - 13/2, 15/2/ 


346221.16 


70.9 


24.03 


41.7(0.1) 


3.6 


0.1) 


7.2(0.2) 


2.6 




H13CO+ 


4-3 


346998.34 


41.6 


60.85 


42.0(0.2) 


1.6 


0.1) 


7.7(0.3) 


2.8 


9 


"CHgOH 


14^ 13A--I4014A- 


347188.28 


254.3 


25.94 


42.7(0.3) 


0.7 


.0.1) 


4.5(0.6) 


0.4 




SiO 


8-7 


347330.63 


75.0 


76.79 


41.0(0.2) 


1.4 


.0.1) 


8.7(0.2) 


2.1 





Notes.(l) Blend with SO at 346528.5. (2) Questionable identification (3) Blend with ^^SO at 337195.0. (4) Blend with ^^SO at 
336590.5. (5) Blend with ^''SO at 351257.2 (6) Blend with CH3CH2CN at 337347.6. (7) Blend with CH3CH2CN at 346947.3. (8) Very 

broad feature with three peaks. (9) Blend with CH3CH2CN at 346983.8. 



The different spatial distribution of molecular gas from 
different species is reminiscent of the chemical differentia- 
tion observed in other hot-cores, like Orion-KL. Most of the 
nitrogen-bearing molecules peak in the Orion hot core while 
most of the oxygen -bearing molecules are found toward 
the c ompact ridge fe. s.. Blake et all Il987l : iBeuther et al.l 
I2OO5I : [Oin ct al. 2010), similar to what observed in G19.61- 
0.23 by Toin et all (|2oTol ). At 7 mm, however, a complex 
oxygen-bearing molecule (acetone) has been detected toward 
the Orion hot core, while two nitrogen-bearing molecules 
(cyanopolyynes) are found in the quie scent cold gas of 
the Orion extended ridge (|Goddi et al.l [2009), similar to 
what observed in G20.08-0.14N. Further higher sensitivity 
and resolution observations especially from the ALMA are 
needed to confirm the peak offsets in G20.08-0.14N. 



3.3 Column Densities and Abundances 

In our observations we detected multiple transitions (> 3) 
from SO2, CH3OH, and HC3N, but owing to spectral blend- 
ing of lines from SO2 and HC3N, we used only CH3OH to 
calculate the rotation temperature and the column density 
by a rotation temperature diagram (RTD). Eight transitions 
of CH3OH have been detected in G20.08-0.14N, containing 
5 ground state and 3 vibrationally excited lines. With the 
assumptions of the local thermodynamic equilibrium (LTE), 
lines being optically thin and gas emission filling the beam, 
the rotation temperature and the beam-averaged column 
density can be determined by l|Goldsmith fc Langeij Il999l : 
iLiu et alll200^: lOin et al.ll2010l) 



. ,N.u Nt . Eu 

ln(_) =ln( — )- — , (1) 

yu ^rot J- rot 
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Figure 3. Sample images from C^*S and SiO in different velocity channels. In each panel, the synthesized beam is shown in the lower- 
right corner. The red " * " and filled triangle indicate the position of three HII region and H2O maser, respectively. Central velocities 
are indicated in each image. Contour levels are all -4, 4, 8, 12, 20, 28, 36, 44 cr .... The oblique dashed line is used to determine the peak 
positional move of core, (a) C'^*S, 1 a noise level is 0.1 Jy beam~^. 



Where Nu is the column density of the upper energy level, 
Qu is the degeneracy factor in the upper energy level, Nt 
is the total beam-averaged column density, Qrot is the rota- 
tional partition function, Eu is the upper level energy in K, 
and Trot is the rotation temperature. By plotting the data 
points from eight transitions of CH3OH according to Equa- 
tion (1) and applying least-square fitting for a straight line, 
a rotation temperature diagram is shown in Figure 6. The 
RTD can be corrected by multiplying Nu/Qu by the opti- 
cal depth correction factor Cr = t/{^ ~ e^), where r is the 
optical depths, which is expressed by (Remijan et al. 2004) 



SKAvTrotgu ' 



(2) 



Where S is the line strength, /j, is the dipole moment, :/(GHz) 
is the rest frequency, k is the Boltzmann constant, and Av 
is the FWHM line width. From the optical depth-corrected 
data, we derived a rotational temperature of 105 ± 29 K 
and the beam-averaged column density of (3.1±2.1)x 10^^ 
cm~^ for the G20.08-0.14 N regio n. 

Following iQin et al] (|2010l ) (their eq.5), the beam- 
averaged column densities of molecules with less than three 
transitions detected could be expressed by 



Ai'T(cm"^) = 2.04 X 10^ 



I (Trot) Qrote J Idv 



for linear molecules {Tiot=Eu)- 



(3) 
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Figure 4. (b) SiO, 1 cr noise level is 0.1 Jy beam ^. 



AfT(cm"^) = 2.04 X 10^ 



(4) 



for symmetric and asymmetric top molecules {Trot=2/3Eu)- 
Where I{Trot) and 1^ are the specific intensities of the spec- 
trum at Trot and of the background continuum, J Idv is the 
integrated intensity of the specific transition in Jy beam~^, 
9a and 6b axe the FWHM beam size in arcsec^. 

In addition, the fractional abundance of a certain 
molecule relative to H2 depends on the column densities 
of the certain molecule and H2 molecule, defined by = 
Nt/N-h^ , where Nt is the beam-averaged column density. 
The optically thin submillimeter dust continuum emission 
has been prov en to be an effective way to determ ine H2 col- 
umn density (|Pierce et al.ll200(]| : iGordonI llOgSl V Assuming 
an average grain radius of 0.1 micron and grain density o f 



the beam-avera ged column density is given by the formulae 
(|Lis et al.lll99ll ') 



ArH2(cm"^) = 8.1 X 10^ 



Q{u)D. ^Jy^^GHz-^ 



(5) 



where T is the mean dust temperature (K), Q(i^) is grain 
emissivity at frequency v, is the peak intensity of the 
continuum, and Q is the beam solid angle. As suming radi- 
ation equilibrium, iGalvan-Madrid et al.l (|2009f ) set the dust 
temperature to 230 K for G20.08-0.1 4N. We adopt Q{u) of 
4 X 10"^ at 340 GHz (|Lis et al.lllQQll ') and dust temperature 
of 230 K for the calculation of the beam-averaged column 
density. The intensity of the continuum peak is 1.33±0.02 Jy 
beam~^. The derived beam-averaged column density of H2 
is (8±0.1) X lO'^^ cm'^, which is reasonably consistent wit h 
those (~10^^ cm^) in NH3 (iGalvan-Madrid et al.l 120091 ') . 



3 g cm-^* and a gas to dust ratio of 100 (|Lis et al.|[l99lh . (3x10^^ cm^) in hot core G327.3-I-0.6 (|Gibb et al.ll2000l) and 
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Figure 5. The integrated intensity maps of eacli species. Contour levels are all 3, 4, 5, 6, 7, 8, 9, 10 cr .... The 1 a noise levels (Jy beam" 
km s^^) of various species are presented in Table 1. 



© 2002 RAS, MNRAS 000,[llfT3l 



8 Jin-Long Xu and Jun-Jie Wang 




100 



200 



300 



400 



500 



E,(K) 



Figure 6. Population temperature diagram of the observed 
CH3OH transitions. The vertical bars mark the l^{N^/g^) er- 
rors from the integrated intensities. The linear least-squares fit 
(solid line) gives a rotation temperature of 105 it 29 K. 



Table 2. The Parameters Derived From Molecular Lines 



Molecule 


Trot 


Nt 




/h2 






(K) 


(cm^ 








SO2 


112 


(2.5 ±0.1) 


X IQi" 


(3.2 ±0.1) X IQ- 


-8 




192 


(3.2 ±0.2) 


X IQis 


(4.2 ±0.2) X IQ- 


-8 


SO 


143 


(3.7 ±0.2) 


X IQi'^ 


(4.8 ±0.2) X 10" 


-7 


C34S 


50 


(5.1 ±0.1) 


X 10" 


(6.5 ±0.2) X 10" 


10 


NS 


71 


(3.3 ±0.1) 


X 10^5 


(4.2 ±0.2) X 10- 


-9 


HC3N 


307 


(7.8 ±0.2) 


X 10^5 


(1.0 ±0.1) X 10" 


-8 


NH2CHO 


100 


(4.4 ± 0.5) 


X 10" 


(5.7 ±0.3) X 10- 


10 




33 


(7.0 ±0.2) 


X 10^'^ 


(9.0 ±0.2) X 10- 


-7 


SiO 


75 


(1.8 ±0.1) 


X lO^s 


(2.3 ±0.2) X 10" 


-9 


H"CO+ 


28 


(2.2 ±0.2) 


X 10" 


(2.9 ±0.2) X 10- 


10 


CH3OH 


105±29 


(3.1 ± 2.1) 


X 10^"^ 


(4.0 ±2.7) X 10" 


-7 


i^CHaOH 


170 


(3.5 ±0.7) 


X 10^5 


(4.5 ±0.8) X 10- 


-9 


HCOOCH3 


165 


(3.5 ±0.3) 


X 10" 


(4.5 ±0.4) X 10" 


-7 



(8.4x10^^ cm^) in hot core G19. 61+0.23 (|Qin et aLlbOlol ). 
Therefore, the derived density of H2 is reliable. The frac- 
tional abundances of the various species relative to H2 are es- 
timated from the beam-averaged column densities as shown 
in the fourth column of Table 2. 



4 DISCUSSION 
4.1 Kinematics 



Recently, 



through 



numerical 



simulation, 

IVazguez-Semadeni et al] (|2009l ) suggest that the for- 
mation of massive stars is coincident with large-scale 
collapse, while low-mass and intermediate-mass stars are 
associated with isolated accretion flows. Previous observa- 
tions with the SMA in hot core molecules (CH3CN, OCS, 
and SO2) and the VLA in NH3 show that G20.08-0.14N 
is surround ed by a smaller scale and a large scale accre- 
tion flows (|Galvan-Madrid et al.|[2009l ). The smaller scale 



accretion flow is across the HIT region A, which may be 
resupplied by the large scale accretion flow. The channel 
maps of SiO and C"^*S in Figures 3-4 exhibit a velocity 
gradient along the SW-NE, suggesting that this is an 
inflow or rotation motions around the HIT region A of 
G20.08-0.14N. 

To explore whether these are two accretion flows in 
G20.08-0.14N, we made the position- velocity (P-V) dia- 
grams (Figure 7) of SiO and C^^S lines across the HIT region 
A with cuts at P.A.=45° and 135°. Since SiO is more eas- 
ily affected by the excitation conditions, the PV diagrams 
of SiO present the complex pattern, which contain several 
velocity components. To analyze in detail each component, 
we divide the whole component into six region. Each has 
been designated alphabetically: Regions D-I. From Figure 7 
(a) and (b), we flnd two velocity components in Region D 
and Region I, which are not seen in the perpendicular di- 
rection P.A. = 135°. The velocity component of Region D is 
from 32.5 km s~^ to 36.0 km s~^, while the velocity com- 
ponent of Region I is from 47.0 km s~^ to 51.7 km s~^. By 
comparin g with previous observations in CH3CN, OCS, SO2 
and NH3 (IGalvan-Madrid et al.]|2009l ). we find that the two 
velocity components of SiO detected in Region D and Re- 
gion I d o not belong to the ve l ocity ranges of the accretion 
flows of iGalvan-Madrid et all (|2009l ). Since SiO is consid- 
ered as a good tracer of outflow (e.g., Schilke et al. 1997b; 
Gueth & Guilloteau 1999; Cesaroni et al. 1999; Beuther et 
al. 2005), we suggest that the two velocity components of 
SiO may arise from an outflow in the NE-SW direction. To 
confirm the interpretation of the two velocity components 
in terms of outfiow, we made the integrated intensity maps 
for each component as shown in Figure 1. In Figure 1, the 
blueshifted and redshifted components are presented as blue 
and red contours. The map clearly shows NE-SW bipolar 
components centered at the peak position of the continuum. 
Hence, we suggest that the velocity gradients seen in Region 
D and Region I are caused by the bipolar outfiow motions. 
The blueshifted and redshifted lobes of the outfiow are in 
a straight line at the NE-SW direction. Klaassen & Wilson 
(2008) sugge sted an outflow with S i O line proflles in G20.08- 
0.14N, while iGalvan-Madrid et all l|2009l ) did not detect the 
outflow in CO lines. Here, we give for the first time the direc- 
tion of the outfiow. A dynamic time scale of the outflow can 
be determined hy td ~ r/v, where v is the maximum flow 
velocity relative to V^ys, and r is the length of the begin-to- 
end flow extension for each lobe. From Figure 1, we derived 
that the lengths of the blueshifted and redshifted lobe S are 
0.15 pc and 0.36 pc, respectively. The obtained average dy- 
namical timescale is about 2.6 x 10'* yr. An H2O is located 
at the redshifted lobe of the outflow, which may be excited 
by the shocks from the outflow. The HII region A is situated 
at the center of the outflow and may drive the coUimated 
outflow. 

Moreover, the velocities of Region F and Region G at 
peaked position are 39.0 km s"^ and 44 km s~\ respec- 
tively, which both may have a velocity gradient of ~ 2.0 
km s~^ with respect to VnyB- The velocity gradient is well 
coincident with the inward ve locity found in the NH3 ab - 
sorption line at smaller scale (|Galvan-Madrid et al.l [20091 ) . 
While the velocity of both Region E and Region H with re- 
spect to Vsys is ~ 4.0 km s~^, which i s well associated with 
the rotation velocity at larger scale iGalvan-Madrid et al.l 
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Figure 7. PV diagrams of SiO and C'''*S. The upper diagrams are for a cut at P.A.=45 (NE-SW), while the below diagrams are at 
P.A.=135 (NW-SE). Cuts are across the position of the HII region A. The black dashed line indicates Vsys- Other blue and red lines are 
used to mark the beginning and end of each velocity component. 



(|2009l ). Additionally, from Figure 7 (c) and (d), we find also 
that the PV diagram of C^^S show a velocity gradient at 
P.A.=135° and P. A -45°. C^'^S is u sed to trace the Kep- 
lerian rotation (|Beuther et al.|[2009l ). (|Galvan-Madrid et al.l 
I2OO9I ) also observed a velocity gradient in molecular lines at 
P. A. =45°, which they interpret as rotation in a torus/disk 
whose plane is oriented NE-SW. Our measurements how- 
ever provide clear indication of an outflow along the same 
NE-SW direction. Therefore, we favor a model where the de- 
tected velocity gradients at P.A.=135° degrees are caused by 
rotation motions (and possibly infall),in the perpendicular 
direction of the outflow. However, in this model we can not 



explain the velocity gradients detected at P. A. =45° degrees 
in the Regions E-H of Figure 7 (a) and Figure 7 (c). 



4.2 Chemistry 

From Table 1, we can see that 11 species were detected. The 
emission peaks for different molecules are located at differ- 
ent positions in the core (Figure 5), hence each molecule 
may be produced via a different mechanism. The individual 
molecules are discussed in the following. 
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4-2.1 Sulfur-bearing Molecules 

The sulphur chemistry is of specific interest because of its 
rapid evolution in warm gas and since the abundances of 
sulphur-bearing species increase significantly with tempera- 
ture, both by ice evaporation and by shock interaction. 

Sulphur dioxide (S02)-Eight transitions of SO2 were 
detected, containing 4 ground-state and 4 vibrationally ex- 
cited lines. S02(19i ig-lSg ^g) at 346.65216 GHz is strongest 
and unblended line. Because SO2 is an inorganic asymmetric 
molecule, adopting a rotation temperature of Trot='2/3Eu, 
we estimated that the column density and abundance of 
SO2 are (2.5 ± 0.1) x lO^'' and (3.2 ± 0.1) x 10"*, respec- 
tively. The derived abundance of SO2 i s consistent with 
those in Sgr B2, G29.26 and G19.61-0.23 jNummelin et al. I 
l200d : iBeuther et a L 2009; ' Oin et al.ll2010l 'l. but larger than 
that in Orion KL ( Beuther et al.l \200^ . It has been re- 
cently suggested that the Orion-KL hot core may be only 
a pre-existi ng density enhancement heated from t he outside 
by shocks JZapata et al.ll201ll : iGoddi et aLll201ll ) . Via the 
shock-induced chemical models, Hartguist et al. (|l980l ) es- 
timated that abundance of SO2 is 3.7 x 10~^^. The highest 
fractional abundance relative to H2 in G20.08-0.14N can be 
explained by grain surface chemistry, not by shock inter- 
action of a molecular outflow with the ambient dense gas. 
Two transitions of the ^*S02 isotopologue were detected. 
In Orion, SO2 dominates the appea rance of the millimeter- 
wave spectrum jSchilke et al.lll997al 'l , accounting for approx- 
imately 28% of all the detected lines. Because of its asym- 
metric geometry, it has a rich spectrum of lines which are 
typically very strong because of the large abundance and 
high dipole moment of the molecule. In G20.08-0.14N, the 
10 transition lines of the detected 41 transition lines be- 
long to SO2 and its isotopologue, which is similar to that in 
Orion. 

Sulphur monoxide (SO)-Sulphur monoxide was de- 
tected in (lO^j^-lOio) transition at 336.55375 GHz. The de- 
rived column density and abun dance is (3 . 7±0. 2) x lO^'' and 
(4.8±0.2) X 10"^ respectively. |Chamie3 (119971 ) has investi- 
gated the formation of SO and SO2. When t^lO'' yr, the evo- 
lution of the sulfur chemistry in the models shows that sul- 
fur monoxide is produced by reaction SH-I-O— >-SO-|-H; When 
tj52xl0'' yr, the reaction S-I-OH— >-SO-|-H is also important 
at 100 K, along with S-I-O2— >-S0-|-0. For increasing core 
temperature, reaction S-|-02^S0-|-0 dominates. Consider- 
ing the SO2 abundance of (3.2±0.1) x 10~*, we obtained the 
SO/SO2 abundant ratio of 15± 0.2 in G20.08-0.14N, which 
can be explained by the evaporated mantles model without 
O2 injection at 200 K. From the model we estimated that the 
age of the G20.08-0.14N hot core is about 8x10^ yr, which 
is less than the dynamical timescale of the outflow identified 
by SiO. In Figure 5, the emission peaks of SO is located to 
the northwest of the continuum peak, which is perpendicular 
to the direction of the outflow. Hence, although the outflow 
can excite the formation of SO from the timescales, SO may 
be produced by reaction SH-f O— >-SO-|-H. After some time, 
we may detect the SO excited by the outflow. 

Carbon monosulfide (CS). The rare carbon-sulfur iso- 
topologue C^^S (7-6) is detected at 337.39646 GHz. The 
estimated column density and fractional abundance of C'^'^S 
is (5.1 ± 0.1) X 10" and (6.5 ± 0.2) x 10"^°, respectively. 
The derived column density of C^^S is similar to the one 



observed in Orion Kl jBeuther etHI |2009|) . [ Beuther et al.l 
(2009) considered that the formation of C'^^S can be suc- 
cessfully explained by gas chemistry models. Moreover, the 
PV diagrams of C'^^S show the inflow and rotation motions, 
further confirming that C^*S should be a better trace of 
rotational disk at very early evolutionary stages. 

4-2.2 Nitrogen-bearing Molecules 

Theoretical models and observations suggest that the 
nitrogen-bearing molecules have higher gas temperatures 
and lower fractiona l abundances on a timescale of^lO^ 
vr lle.g.. Blake et al.l [l987l: iRodgers fc CharnlevI 12001.20031 : 



iMiao et al.lll995l : lKuan fc Snv der"l996': 'Oin et al."20ldl). In 
the G20.08-0.14N hot core. [Calvan-Madrid et al., (20091 ) de- 
tected NH3 and CH3CN molecules. They derived that the 
fractional abundances of both molecules are about 5xl0~^ 
and 5x 10"'' to 2x 10"^. NH3 and CH3CN are the best trac- 
ers of hot cores (le.g.. Wilner et al. l ll994l : IWright et a"Llll996l : 
I Wilson et al.ll2000[). NH3 may ori gin from the evaporation of 
grain mantles l|Pauls et al.lll983l ). which drives a nitrogen- 
rich chemis try and produces mu ch complex nitrogen-bearing 
molecules (|Caselli et al.l Il993l ). Here, we detected three 
nitrogen-bearing molecules, including cyanoacetylene, nitro- 
gen sulfide, and formamide toward the G20.08-0.14N hot 
core. 

Cyanoacetylene (HC3N). We detected the two ground- 
state and three vibrationally excited state lines of HC3N. 
HC3N (38-37e) at 346.45573 GHz is strongest and unblend. 
Using this transition line, we obtained the HC3N column 
density and fractional abundance of (7.8 ± 0.2) x 10^^ and 
(l.OiO.l) X 10"*, respectively. The derived fractional abun- 
dance of HC3N is one order of magnitude larger than that in 
Orion KL and Sgr B2 l|Beuther et al.ll2009l : lNummelin et al. I 
ioo^). The fractional abundance of HC3N with respect to 
the NH3 is 0.02. 

Nitrogen sulfide (NS) Nitrogen sulfide at 346.22116 GHz 
is the first detected in the J — 15/2 — 13/2 transition. NS 
is a relatively simple species, which provides a good test 
for coupled chemistry models of nitrogen and sulfur. The 
derived column density of NS is (3.3 ±0.1) x 10^^. The frac- 
tional abundance of NS relative to molecular hydrogen is 
(4.2 ± 0.2) X 10"", which is larger than about 10"^^ found 
by some recent g as-phase chemis try models dev eloped for 
quiescent clouds (|Lee et al.l 1 19961 ). In addition. IViti et al.l 
(2001) predicted that the abundance of NS is significantly 
enhanced by shock. NS does appear to follow the distri- 
butio n of SO and CH3OH in TMC-1 (|McGonagle fc Irvine! 
I1997I ). In G20.08-0.14N, the emission peak of NS is similar 
to the distribution of HC3N and SiO. The emission of the 
SiO shows the outflow. Hence, we suggest that the formation 
of NS may be related to the outflow of G20.08-0.14N. 

Formamide (NH2CHO) We also detected NH2CHO at 
336.13688 GHz. Assuming a rotation temperature of Trot 
= 2/3Eu, the estimated column density of NH2CHO is 
(4.4 zfc 0.5) X 10^"^, which is close to val ues in Orion K L 
(jSchilkeet al.lll997al ) and in G19.61-0.23 (|Qin et al.ll20i(il ). 
Because the abundance of NH2CH O is 10 times less than 
that of the HNCO in G19.61-0.23, lOin et al.l (|2010l ) sug- 
gest that NH2CHO can be only formed by successive hydro- 
genati on of the HNCO on grain surface. iTielens fc CharnlevI 
(|l997h considered that NH2CHO is most likely formed 
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by atom addition to HCO on grain mantl es and then 
evaporated. Moreover, iBernstein et alj (| 19951 ) found that 
NH2CHO is formed through UV photolysis and heating of 
a H2O:CH3OH:CO:NH3 = 100:50:10:10 mixtu re. If NH2CHO 
is formed via UV emission, iGibb et al.l l|2000l ) suggested that 
the abundance of CH3OH is higher enough than that of NH3. 
For G20.08-0.14N, the peak position of NH2CHO with re- 
spect to the continuum is similar to the H^'^^CO"'', and the 
abundance ratio of CH3OH is lower than that of NH3, in- 
dicating that the formation of NH2CHO m ay be similar to 
the suggestion of lTielen7 fc CharnlevI (|l997l V 



4-2.3 Oxygen-bearing Molecules 

Comparing with the nitrogen-bearing molecules, the oxygen- 
bearing molecules have lower gas temperatures and higher 
fractional abundances on a timescale of^lO* yr (Blake et 
al. 1987; Rodgers & Charnley 2001,2003; Qin et al. 2010). 
CH3OH and HCOOCH3 were detected in G20.08-0.14N, 
which are the typical oxygen-bearing molecules. CH3OH is 
considered as the precursor of the oxygen-bearing molecules, 
which drives an oxygen-rich chemistry and produce new 
molecules such as HCOOCH3. 

Methanol (CH3OH). Methanol has been detected in 
CH3OH and its isotopologue ^^'CHgOH. There are eight 
transitions of the CH3OH with upper level energies of 61- 
489 K including 5 ground state and 3 vibrationally ex- 
cited lines. Through the RTD fit, the derived rotation tem- 
perature of the CH3OH is larger than 100 K, suggest- 
ing that they originated from warm gas environments. In 
addition, the column density and fractional abundance of 
CH3OH are (3.1 ± 2.1) x lO" and (4.0 ± 2.7) x 10"^ re- 
spectively. The obtained fractional abundance of CII3OH is 
one ord er of magnitude lower than that in Orion KL and 
Sgr B2 (jBeuther et al.ll2009l : iNummelin et al. II2OO0I), but is 

and 
and 



consistent with tha t in G19.61-0.2 3 (iQin et al.l 



2OIOI') 
20101 ) 



IRAS 20 1 26-K41 04 (|Xu et all [20]! '). lOin et al.1 

IXu et al.1 \2Qli ) suggest that the higher fractional abun- 
dance of CH3OII with respect to II2 can be explained by 
grain surface chemistry. The grain surface chemistry mod- 
els pr edict that the abundan ce of CH3OH is larger than 
10~^ (|van der Tak et al. I2OO0I ), while t he gas chemical m od- 
els predict an abundance of< 10~® (|Lee et al.l 1 19961 ). In 
G20.08-0.14N, the higher fractional abundance of CH3OH 
cannot be interpreted simply by gas-phase chemical reac- 
tions, and may be originated from grain surface chemistry. 
Th e ^^C/"C ra t io fro m the CH3OH and "CH3OH is ~ 
89. iMilam et al] (|2005l ) concluded that the ^^C/^^C ratio 
obtained from CO indicates a gradient with Galactic dis- 
tance of ^^C/^^C=5.41D-|-19.03, where D is distance from 
the Galactic center in kpc. Using the ^^C/^^C ratio from 
the CII3OH, we derived that the distance of G2 0.08-0.14N 
is 12.9 kpc . which is close to valu es reported bv lFish et al.l 
l|2003t ) aud i Anderson et al.l (|2009l ). 

Silicon monoxide (SiO). SiO J=8-7(v=0) was detected 
at rest frequencies of 347.33063 GHz with upper level en- 
ergy of 75 K. The column density and fractional abundance 
derived from a single line of SiO are (1.8 ± 0.1) X 10^^ and 
(2.3 ± 0.2) X 10~®, respectively. The derived column den- 
sity of SiO is an order of magnitude higher than that of 
Orion KL and G29.96 (|Beuthcr ct al, ,2005 ). One explana- 
tion for the SiO abundance enhancements in shocked gas is 



that SiO is produced by destruction of g rain cores in shocks 
(C aselli et al.1 1 19971 : ISchilke et al.lll997bl ). Another possibil- 
ity is that SiO is embedded in icy grain mantles, which are 
evaporated. Hence, we suggest that SiO may be produced 
by destruction of grain cores in the shocks of the outflow. 

Formyl ion (HCO+) and Methyl formate (HCOOCH3). 
We detected the (4-3) transitions of the H^'^CO"'" isotopo- 
logue. The estimated column density and fractional abun- 
dance ofHi3CO+ are (2.2±0.2) x lO" and (2.9±0.2) x 10"^°. 
In Figure 5, the emission of HCO^ shows an elongated struc- 
ture with two cores. Since the HCO^ is blended with the line 
from CH3CH2CN at 346983.8 MHz, the latter could be re- 
sponsible for the complex structure of HCO^. HCOOCH3 
is a heavy asymmetric rotor with hindered internal rotation 
of the methyl group. Two transitions of HCOOCH3 were 
detected with the relatively weak emission. Using the tran- 
sition of HCOOCH3 at 347.47824 GHz, we obtained the col- 
umn density and fractional abundance of (3.5 ± 0.3) x 10^^ 
and (4.5 ± 0.4) x 10"^, respectively. Our obtained column 
density of HCO OCH3 is consistent with that of Sgr B2(N) 
(|Liu et al.|[200ll ). but larger than that (~10~*) produced 
by gas phase chemistry. 

In the Orion KL hot core, because the hot molecular 
gas is not associated with any self- luminous mi llimeter, ra- 
dio, or embe d ded i nfrared source, iGoddi et al.1 (,201 1 ) and 
IZapata et al.l (|201ll ) suggested that the Orion-KL hot core 
may be only a pre-existing density enhancement heated from 
the outside by shocks from a molecular outflow. On the other 
hand, in G20.08-0.14N the molecular gas is associated with 
a millimeter source and an HII region and clearly shows star 
formation activity. This indicates that the G20.08-0.14N hot 
core is heated by a protostar forming at its center. Based on 
the chemical models discussed here, we concluded that the 
age of the G20.08-0.14N hot core should be about 10*-10^ 
years. 



5 SUMMARY 

The Submillimeter Array observations toward the high-mass 
star- forming region G20.08-0.14N are presented in the sub- 
millimeter continuum and in the molecular line transitions. 
From the SMA data, 0.9 mm continuum emission reveals an 
extended structure which is associated with an H2O maser 
and three UC HII and HC HII regions. Forty-one molecu- 
lar transitions were detected related to 11 molecular species, 
including SO2, SO, C^^S, NS, C'^'O, SiO, CH3OH, HC3N, 
H^='CO+, HCOOCH3, NH2CHO and their isotopic species. 
In addition, six molecular transitions are unidentified. 10 
transition lines of the detected 41 transition lines belong to 
SO2, which dominates the appearance of the submillimeter- 
wave spectrum. The channel maps of C'^'^S and SiO show 
velocity gradients. In their PV diagram, C^^S emission show 
rotation motions, while SiO not only present two rotation 
motions in smaller and larger scales respectively, but also re- 
veal for the first time a coUimated outflow along the NE-SW 
direction. The average dynamical timescale of the outflow is 
about 2.6 X 10* yr. An HII region is situated in the central 
position of the outflow, which may drive the collimated out- 
flow. Eight transitions of CH3OH are unblended, we derived 
the rotational temperature and the column density of 105 K 
and 3.1 x 10^^ cm^^ for CH3OH lines, respectively, further 



© 2002 RAS, MNRAS OOO.fTHTSl 



12 Jin-Long Xu and Jun-Jie Wang 



indicating that a hot core coincides with G20.08-0.14N. The 
hot core is heated by a protostar at its center with the age 
of about 10''-10^ The emission peaks of different molecules 
are located at different positions of the hot core. By compar- 
ing the abundances of different species with chemical models 
and previous observations of other hot cores, we concluded 
that each species may be produced with a different mecha- 
nism. Nitrogen sulfide (NS) is for the first time detected in 
G20.08-0.14N. 
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